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During epitaxial crystal growth a pattern that has initially been imprinted on a surface approxi- 
mately reproduces itself after the deposition of an integer number of monolayers. Computer simu- 
lations of the one-dimensional case show that the quality of reproduction decays exponentially with 
a characteristic time which is linear in the activation energy of surface diffusion. We argue that 
this life time of a pattern is optimized, if the characteristic feature size of the pattern is larger than 
{D / F) 1 ^ d+2 \ where D is the surface diffusion constant, F the deposition rate and d the surface 
dimension. 



PACS numbers: 68.55.-a, 05.50.+q, 81.15 

Modern techniques of manipulating crystal surfaces al- 
low to imprint structures down to the atomic size on 
them. With the tip of a tunneling microscope one can 
arrange adsorbed atoms in a pattern. In heteroepitaxial 
growth a two dimensional array of quantum dots may 
form on the nanometer scale. With masking techniques 
arbitrary patterns with features as small as a microm- 
eter along the surface and atomic size perpendicular to 
it can be fabricated. If one buries such a pattern un- 
der an overlayer growing in layer-by-layer mode, the pat- 
tern will approximately reproduce itself periodically at 
the completion of each layer. The question arises, how 
this propagation of a pattern is influenced by the growth 
conditions. A theoretical understanding of such tempo- 
ral correlations began to emerge only recently [0-|5|. In 
the following the propagation of a pattern will be dis- 
cussed for the simplest case that the influence of elastic 
strain on surface diffusion is negligible and that diffusion 
across a step edge from an upper to a lower terrace is not 
inhibited by Ehrlich-Schwoebel barriers ||. 

Experimentally, it has long been known that layer-by- 
layer (or Frank-van der Merwe) growth persists up to a 
time t, after which the characteristic growth oscillations 
are damped out and the surface becomes rough. Re- 
cently, it has been shown that t depends on the growth 
conditions i.e. the surface diffusion constant D 

and the deposition rate F, with a power law (D/F) s . 
Clearly, a pattern can at most survive as long as the sur- 
face grows layerwise. In fact the life time of a pattern is 
much shorter than t, as we are going to show: It depends 
on D/F only logarithmically rather than with a power 
law. 

The life time of a pattern in layer-by-layer growth de- 
pends also on the feature size r. We shall identify the 
length scale that is important in the context of pattern 
decay as £ ~ {D / F) 1 / ( - d+2 \ For feature size r > l , the 
life time increases rather abruptly. 



Propagation probability. To simplify the discussion, we 
assume that the grown crystal is simple cubic [jioj] and 
that overhangs and defects can be neglected. Then, the 
surface configuration at a given time t can be represented 
by a height function h(x,t) with x indexing the lattice 
sites on a (i-dimensional substrate. In the following, time 
will be measured in numbers of deposited monolayers. 

When the growth process is initiated at time t = with 
some structure h(x, 0), consisting of e.g. islands of a par- 
ticular shape on an otherwise flat substrate (a pattern ) , 
new islands will form on top of the initial ones, and all 
of them will expand laterally due to the attachment of 
atoms at their edges. Therefore the pattern will be de- 
formed during this early stage of growth. However, the 
original pattern will nearly be reproduced after the depo- 
sition of one monolayer, because the new layer nucleates 
preferentially near the centers of islands in the previous 
layer. As this correlation extends over long times [jlj, one 
expects the approximate reproduction of the pattern at 
later times as well. It makes sense to ask, which fraction 
of the pattern is propagated through t monolayers. We 
call this fraction the propagation probability p(t). It is 
given by 



P(t) = (H £ft(x,s),Mx,0)+s) 
s=l 



(1) 



The brackets denote averaging over different lattice sites 
x and different realisations of the growth process. Sij = 1 
if i = j and otherwise denotes the Kronecker delta. 

By defining p(t) in this way, we measure the deter- 
ministic reproduction of the initial pattern h(x,0), seen 
stroboscopically after deposition of 1, 2, . . . monolayers in 
the comoving frame. It is deterministic in the sense that 
a given site x is counted as 'surviving' after time t only 
if it has survived through all previous times 1, . . . , t. Of 
course, the height might also regain its initial value by 
chance once it has left it, but this is a stochastic process 
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and contributes to the noisiness of the pattern. Thus 
the survival of information is not described properly by 
the propagation probability (Q), but by an appropriate 
entropy. It will be dealt with in a longer paper. 

If the probability of return to the initial height (in the 
comoving frame) at time t was considered, irrespective of 
the values at intermediate times, one would probe a two- 
point function. Such a two-point function is expected 
to decay algebraically for long times. More precisely, it 
should scale like the value of the height distribution func- 
tion at the average height. For self-affine surfaces the 
width of the height distribution increases like , Thus 
the probability to recover the initial height after the de- 
position of t monolayers, should decay as t~@. 

In contrast to this, the t— point function p(t) decays 
exponentially, see fig. [j]. Although the exact evaluation 
is non-trivial, this result is easily made plausible: The 
fraction p(t+l) of surviving sites at time t+1 equals the 
number of surviving sites at time t, p(t), times the prob- 
ability to propagate to the next layer. Assuming that 
this probability can be identified with p(l), independent 
of the surface configuration, which at time t of course 
differs from the initial one, the exponential decay 

p(t) = exp(-i/i c ), (2) 

follows immediately. We shall show below that the prop- 
agation probability does depend on the surface configura- 
tion. However, this dependence is so weak, that the sur- 
face evolution during the life time t c of a pattern hardly 
affects its exponential decay (^|). 

The main purpose of this Letter is to investigate the 
dependence of t c on the microscopic growth parameter 
D/F and the feature size of the pattern. Let us first 
discuss two limit cases. For D/F — > 0, the sites are 
not coupled by diffusion and the appropriate description 
of the gro wth process is given by the random deposition 
model [[llj , where atoms are deposited randomly onto the 
substrate and remain at the deposition site forever. In 
this model, the fraction Sj (t) of surface sites in layer j af- 
ter time t is a Poisson distribution, Sj(t) — exp(—t)P/j\. 
Therefore, p(l) = Si(l) = exp(-l), so that t c = 1. 

For D/F ^> 1, a monotonous increase of t c as function 
of D/F is expected. When the initial 'pattern' is simply a 
flat surface, one expects t c — > oo for D/F — > oo, because 
layer-by-layer growth persists forever for infinitely high 
diffusion constant. The computer simulations of molecu- 
lar beam epitaxy (MBE) on a one-dimensional substrate, 
to be presented in the next section, confirm this picture 
and show a dependence 

t c ~ log(£>/F) (3) 

with a cutoff at t c = 1 for small D/F. 

In the following we shall discuss three different initial 
patterns: (1) A completely flat surface, (2) a rough sur- 
face as it evolves from the flat one after time t, when 



the oscillations due to layer-by-layer growth have died 
out, and (3) a periodic modulation of the surface with a 
fixed feature size. The first arises as a natural limit of 
a pattern with a characteristic feature size r — > oo. The 
second represents the simplest generic configuration for 
which the growth kinetics has no periodic time depen- 
dence any more. Both will be used to study the pattern 
decay process systematically in the next section. Finally, 
the third choice will lead to the optimization condition 
for pattern survival and will be studied afterwards. 

Model and simulation results. Atoms are deposited 
onto a one-dimensional substrate of typical size L = 10 4 
with a rate of F atoms per unit time and area. Atoms 
with no lateral neighbour are allowed to diffuse with dif- 
fusion constant D. Atoms with lateral neighbours are 
assumed to be immobile, so that e.g. dimers are immo- 
bile and stable. Growth commences with a flat substrate, 
h(x, 0) = for all sites x = 1, . . . , L. (The other initial 
configurations will be discussed below.) On deposition 
at x, h(x, t) is increased by one. 




FIG. 1. Propagation probability p(t) of a flat surface as 
function of number of deposited monolayers t for different 
values of D/F = 0, 10°, 10 1 , . . . , 10 s from left to right. The 
characteristic time of the exponential decay increases with 
increasing D/F. 



Fig. [j] confirms the exponential decay of the propa- 
gation probability. The initial configuration survives the 
better the higher the value of D/F, as expected. The 
life times t c (D/F) are shown in fig. ||. For small D/F 
the life time approaches 1 as derived above for random 
deposition, and for large D/F it increases like \og(D/F). 

The results for a rough surface as initial pattern are 
very similar. The propagation probability (^) was aver- 
aged over 100 different initial patterns, all obtained by 
depositing 50, respectively 200 monolayers on a flat sub- 
strate. For the values of D/F considered here the peri- 
odic oscillations of the surface morphology have stopped 
by then, as shown in ||. Interestingly, the life time of a 
rough surface pattern is shorter than that of the flat sur- 
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face, see fig. but we could not observe any difference 
between the rougher surface (200 monolayers) and the 
less rough one (50 monolayers). This will become plau- 
sible, when considering, how the life time of a pattern is 
affected by the feature size. 
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FIG. 2. Pattern life time t c as function of D/F for a 
flat (o) and a rough (□) surface. The solid line indicates the 
minimum life time t c = 1 for D/F — » 0. The dotted lines 
represent logarithmic fits to the last five decades of D/F. 

Periodic patterns. In this section, we study the depen- 
dence of the life time on the typical feature size of an 
artificially prepared initial configuration. To this end, it 
is useful to recall the different length scales associated 
with ideal MBE: 

The island distance or diffusion length £ is a function 
of D/F: 



£~(D/F)- 



(4) 



The exponent 7 depends on the substrate dimension, on 
the size of the critical nucleus and the possible fractality 
of the islands |^,^,^2| . Its numerical value for the simu- 
lations presented here is 7 = 1/4. 

The only dimensionless length scale £ which can be 
constructed from the dimensionful parameters D and F, 



(D/F) 



l/(d+2) 



(5) 



Physically, this length scale comes from comparing the 
diffusion time to the adatom arrival time on an area of 
size l d p|,[i3||. £q and i are submonolayer quantities. 

Finally, the layer coherence length £ depends on D/F 
like 



l~ (D/F) 4 i/( 4 - d \ 



(6) 



see . This is not a submonolayer quantity, as £ appears 
as the typical length scale after the oscillation damping 
time, i.e. after deposition of t monolayers. 



In order to study the length scale dependence of the 
life time, we use a periodically modulated surface 



h(x,t = 0) = 6 (sin(7ra/r)) 



(7) 



where &(x) = 1 if x > and otherwise is the Heavyside 
function. 
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FIG. 3. Life time t c as function of the feature size r of 
the initial surface modulation for D/F = 10 4 , 10 5 , 10 6 , 10 7 . 
Around a characteristic feature size depending on D/F, the 
life time increases rapidly. 



The measurements of the life time as a function of the 
initial wavelength r for different D/F, see fig. |^, show 
that for r greater than a characteristic value depending 
on D/F, survival is strongly enhanced. Above this value, 
the life time depends only little on the feature size. For 
small r, the life time seems to saturate for increasing 
D/F. These findings can be understood in the following 
way. 

The mechanism of transporting the memory of the sur- 
face structure from one monolayer to the next, is that 
nucleations take place near the center of islands that 
have already formed one layer below. If the feature 
size r is chosen so small that nucleations cannot take 
place on top of the pattern, this mechanism is suppressed 
and consequently the life time of the pattern is reduced. 
To suppress nucleations, the distance between sinks for 
adatoms, i.e. the feature size r, has to be chosen so small 
that a freshly deposited adatom diffuses to the nearest 
sink (i.e. a distance r), before the next atom is deposited 
within the area ~ r d . This is the case for r < £q. Hence, 
£0 should be the length scale found in fig. [| 

The scaling plot fig. [| shows that the characteristic 
length above which the survival is prolonged, scales like 
(D/F)°- 32±0 01 , which suggests a length scale propor- 
tional to (D/F) 1 / 3 . This is in accordance with the argu- 
ment given above for the characteristic length scale being 
£0 . To assure that the characteristic length scale in fig. || 
is not £, which for the parameters studied here also is pro- 
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portional to (D/F) 1 / 3 , we studied epitaxial growth with 
a critical nucleus of 2 instead of 1, which influences £, but 
not £q. This analysis shows that indeed £q is the charac- 
teristic length scale. Details on this will be published in 
a longer paper. 
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FIG. 4. Data from fig. ^| with normalized life times, and 
feature size r scaled with (D/F) - 32 . 

The saturation of t c for small r as function of D/F 
can now easily be understood: If r <C £q, the memory 
of the initial periodic pattern is destroyed already after 
the deposition of one monolayer. The surface will be half 
filled with islands which have a typical size £ much larger 
than r. Then the fraction of sites which propagated by 
just one lattice constant is about 50%, irrespective, how 
large D/F is. 

The faster decay of a rough compared to a flat surface 
as initial pattern can be made plausible with the follow- 
ing reasoning: The feature size of rough surfaces may be 
identified with the typical terrace size I. For the simula- 
tions presented here, where the critical nucleus was 1, 1 is 
smaller than £q. Therefore the faster decay of the pattern 
is consistent with our findings for periodic patterns. 

Conclusions and outlook. In conclusion we have shown 
that a pattern decays exponentially fast with a life time 
proportional to log(D/F). With the Arrhenius law, D ~ 
cxp(— E/ksT), the life time decreases linearly with the 
energy barrier E for surface diffusion. The life time of 
a pattern is optimal if the feature size of the pattern is 
larger than (D/F) 1 '^ 2 ^. 

An important extension of this work would be the 
study of the two-dimensional case. Whereas it is nat- 
ural to expect an exponential decay of the propagation 
probability, the dependence of the life time on the micro- 
scopic growth parameters is an open question. 

In this paper, we neglected barriers for interlayer 
transport (Ehrlich-Schwoebel barriers) The memory 
mechanism will be enhanced by them, but the instabil- 
ity JlJ] associated with them will tend to make pattern 



reproduction worse. The competition between these two 
effects is well worth studying as in many materials inter- 
layer transport in inhibited. 
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